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Estrogen receptor a polymorphism in boys with 
constitutional delay of growth and puberty 

Purpose: There were a lot of reports regarding associations of polymorphisms in the 
estrogen receptor a (ESR1). with many disorders. But, those with constitutional delay 
of growth and puberty (CDGP) are not known. Our aim is to find out any association 
between CDGP and ESR1 . 

Methods: In a total of 27 subjects, we compared 7 CDGP patients with 20 healthy 
controls with their heights and sexual maturity rates were within normal range. We 
selected three single nucleotide polymorphisms from intron 1 of ESR1 (rs3778609, 
rsl 2665044, and rs827421 ) as candidates, respectively. 

Results: In genotype analyses, the frequency of G/G genotype at rs827421 in intron 1 
of ESR1 was increased in CDGP boys (P=0.03). 

Conclusion: The genetic variation of ESR1 can be a contributing factor of tempo of 
growth and puberty. 

Keywords: Eestrogen receptor alpha, polymorphism, single nucleotide, Puberty 



apem 



Byung Ho Kang, MD, 
SoYoun Kim, MD, 
Mun SukPark, MS, 
Kyung Lim Yoon, MD 
KyeShikShim, MD, PhD 

Department of Pediatrics, 
Kyung Hee University School of 
Medicine, Seoul, Korea 



Introduction 
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Growth and puberty are influenced by genetic and environmental factors. The heritability 
of growth and puberty is reported about 75-90 and 50-80%, respectively. Even though the 
large proportion of them is affected by genetic factors, but it is difficult to elucidate a major one 
because of complex polygenic traits . 

The boys with constitutional delay of growth and puberty (CDGP) usually have short stature 
due to slow tempos of growth and are referred to pediatric endocrinologist. But, it may be a 
challenge to differentiate CDGP from pathologic short stature before onset of puberty, and 
genetic factors involved in CDGP are not clarified yet 7 9) . 

There were several reports about associations of estrogen receptor a (ESR1) polymorphisms 
with bone mineral density, breast cancer, height, and age of menarche 1 . Also, the abundant 
amount of ESR1 is present in epiphyseal plates of long bones 19 '. Therefore, we hypothesized 
that ESR1 can be a factor that controls the tempo of growth and puberty. 

The purpose of this study is to determine whether polymorphisms in ESR1 are associated 
with CDGP in boys. 
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Materials and methods 

1. Subjects 

From July 2008 through August 2012, we recruited 27 boys who visited Kyung Hee 
University Hospital at Gangdong for growth evaluation. Seven patients diagnosed with 
CDGP, which is defined as a condition in which the height of individual is at or below the 
3rd percentile in the same sex and age group according to the Korean population-based 
reference 20 ', annual growth velocity is greater than 4 cm, with delayed puberty and bone age, 
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but, without evidence of systemic, endocrine, nutritional or 
chromosomal abnormalities 7 '. To differentiate CDGP from 
hypogonadotropic hypogonadism or hypopituitarism, they 
underwent combined pituitary function tests 8 ' 9 '. Twenty boys 
with normal height and pubertal range were grouped as a 
control. 

Height (cm) was measured using standardized equipment 
(Harpenden stadiometer Ltd., Crymych, UK). The height 
standard deviation score was calculated by (measured height - 
mean height)/standard deviation. The bone age was assessed by 
radiographs of the hand and wrist, using a method of Greulich 
and Pyle 21 '. Midparental target height was calculated by (Father's 
height [cm] + Mother's height [cm] + 13)/2 in male. 

All subjects and parents provided informed consent before 
the study commencement. This study was approved by the 
clinical research ethics committee at Kyung Hee University 
Institutional Review Board (KHNMC- IRB 2008-019). 

2. Genotyping 

Genomic DNA was isolated from peripheral blood samples 
of subjects using DNA extraction kit (Qiagen, Venlo, The 
Netherlands) according to the manufacturer's recommendation. 
The genomic DNA was amplified with polymerase chain 
reaction (PCR) with forward and reverse primer pairs 
and standard PCR reagents in 10 microliter reaction 
volume, containing 10 ng of genomic DNA, 0.5 pM of each 
oligonucleotide primer, 1 microliter of 10X PCR buffer, 250 
uM dNTP (2.5 mM each) and 0.25 unit Taq DNA Polymerase 
(5 unit/uL) (iNtRON Biotechnology, Seongnam, Korea). The 
PCR reactions were carried out as follows: 5 minutes at 95 °C 
for 1 cycle, and 35 cycles at 95 °C for 30 seconds, 55 °C for 30 
seconds, 72 °C for 30 seconds, followed by 1 cycle of 72 °C for 
10 minutes. After amplification, the PCR products were treated 
with 1 unit each of shrimp alkaline phosphatase (SAP) (USB 
Co., Cleveland, OH, USA) and exonuclease I (USB Co.) at 37 °C 
for 75 minutes and 72 °C for 15 minutes to purify the amplified 
products. One microliter of the purified amplification products 
were added to a ready reaction mixture containing 0.15 pmol 
of genotyping primer for primer extension reaction. The primer 
extension reaction was carried out for 25 cycles of 96 °C for 
10 seconds, 50 °C for 5 seconds, and 60 °C for 30 seconds. The 
reaction products were treated with 1 unit of SAP at 37 °C for 1 
hour and 72 °C for 15 minutes to remove the excess fluorescent 
dye terminators. One microliter of the final reaction samples 
containing the extension products were added to 9 microliter of 
Hi-Di formamide (ABI, Foster City, CA, USA). The mixture was 
incubated at 95 °C for 5 minutes, followed by 5 minutes on ice 
and then analyzed by electrophoresis in ABI Prism 3730x7 DNA 
analyzer. Analysis was carried out using Genemapper ver. 4.0 
(Applied Biosystems, Foster City, CA, USA). Table 1 shows the 
primer sets and Tin used for the reactions. 



3. Statistical analysis 

Statistical analysis was performed using SAS 8.02 (SAS 
Institute Inc., Cary, NC, USA). The demographic characteristics 
of the control and patient group were compared with Kruskall- 
Wallis test. Each result was presented as the mean±standard 
deviation. The value of P < 0.05 is considered as statistically 
significant. The frequencies of allele and genotype, and the 
departures of the genotype distribution from Hardy- Weinberg 
equilibrium for each single nucleotide polymorphism (SNP) 
were analyzed using the chi-square test or Fisher exact test. 
Linkage disequilibrium was calculated with the Haploview 
ver. 3.2 (Broad Institute, Cambridge, MA, USA). The genotype- 
specific risks were estimated as the odds ratios with associated 
95% confidence intervals using conditional logistic regression 
analysis. 

Results 

1. Subjects 

The chronologic age of boys with CDGP is 14.3±0.5 years and 
controls is 14.3±0.3 years. The bone age (12.0±0.5 years) was 
less than chronological age in CDGP patients with significant 
statistical difference (P < 0.05). However, there were no 
significant differences between bone age (14.0±1.0 years) and 
chronological age in the control group. Mean height standard 
deviation score (SDS) in CDGP group was less than that of 
the control group (P< 0.05), but there were no differences in 
midparental target height SDSs between CDGP and control 
group. For sexual maturity rates, data showed Tanner stage 1 in 
CDGP patients and 3 in control group (Table 2). The pituitary 
function tests were performed in all CDGP patients, and the 
results were within normal limits (The data was not shown). 

2. Distribution of ESR1 polymorphisms 

Three SNPs were analyzed; rs3778609, rs 12665044 and 
rs827421 in intron 1 of ESR1. The rs3778609 genotypes were 
C/C, 42.86%; C/T, 28.57%; and T/T, 28.57% in CDGP patients, 
and 50.0%, 45.0% and 5.0% in controls. The rsl2665044 

Table 1. Primer sets and Tm for polymerase chain reaction and 
sequencing reactions 

ESR1 (rs3778609): 55°C 

Forward: ATGAGCATATTTTC 

Reverse: CTGTCCACATTCTTTGACG 
ESR1 (rsl 2665044): 55 °C 

Forward: ATGGCAGATTTGGGAATGTC 

Reverse: ATGTAAAGGGTCGGATACAA 
ESR1 (rs827421):55°C 

Forward: TTTTTTAACTTGAGTTTCCCGG 

Reverse: TACACCTCATAACAACTAGE 

ESR1 , estrogen receptor a; rs, reference sequence. 
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genotypes were C/C, 42.86%; C/T, 28.57%; and T/T, 28.57% 
in CDGP patients, and 45.0%, 45.0%, and 10.0% in controls. 
The rs827421 genotypes were A/ A, 14.29%; G/A, 42.86%; and; 
G/G, 42.86% in CDGP patients, and 35.0%, 60.0%, and 5.0% 
in controls. Therefore, the frequency of G/G genotype at the 
rs82742 1 in intron 1 of ESR1 was increased in CDGP patients 
with odds ratio of 19.15 (/"-value permutation=0.03) (Table 3). 
In haplotype frequencies, there were no significant differences 
between cases and controls (Table 4). 

Discussion 

CDGP is a condition with short stature due to slow tempo 
of growth and puberty without any evidence of systemic, 
endocrine, nutritional, or chromosomal abnormalities. 
Although the heights of children with CDGP are lower than 
3rd percentile before onset of puberty, their adult heights can be 
in the range of target height. Therefore, in those children, many 



Table 2. Demographic characteristics of subjects 



Characteristic 


Case (n=7) 


Control (n=20) 


P-value 


Age (yr) 


14.25+0.50 


14.33+0.33 


NS 


Bone age (yr) 


12.00+0.50 


14.00+1.00 


<0.05 


Height (cm) 


145.50+2.30 


1 56.70±3.80 


<0.05 


Height SDS 


-2.02+0.38 


-0.72±0.63 


<0.05 


MPH (cm) 


1 70.40±4.80 


172.30+2.70 


NS 


MPH SDS 


-0.52±0.80 


-0.18+0.45 


NS 


Testicular size (mL) 


3.00±0.50 


6.00+1 .00 


<0.05 


Tanner stage of pubic hair 


1 


3 


<0.05 



Values are presented as mean±standard deviation. 

SDS, standard deviation score; MPH, midparental target height; NS, not 

significant. 



interventions and modalities of promoting growth are usually 
unnecessary 1,3 '. 

To differentiate CDGP from other pathologic short stature, 
for example, subcategory of idiopathic short stature or 
hypogonadotropic hypogonadism, many studies focused on 
the genes and hormonal evaluations related with pubertal 
timing were done 2 ' 4 " 7 '. But, it is a challenge because of similar 
phenotypes, overlapping of hormonal levels, and still 
unidentified genetic factors of them 7 '. 

In aspect of phenotype, the complete discrimination between 
CDGP and hypogonadotropic hypogonadism may be possible 
after puberty . In hormonal studies, basal serum levels of 
luteinizing hormone, follicle stimulating hormone or inhibin B, 
stimulation tests with gonadotropin-releasing hormone (GnRH) 
or agonist and human chorionic gonadotropin were reported to 
be useful . 

The genetic studies about pubertal timing or CDGP are more 
complicated 1 ' 4 '. Although the inheritance is likely complex, 
several studies to find out major genes related with them were 
committed because some predisposing genetic factors have 



Table 4. Estimated haplotype frequencies of intron 1 of ESR1 in 
CDGP cases and controls 



Haplotype* 


Case (%) 


Control (%) 


OR (95% CI) 


/Tor 
interaction 


CCA 


5(35.71) 


26 (65.0) 


0.28 (0.08-1.04) 


0.06 


CCG 


3(21.43) 


2 (5.0) 


4.78 (0.68-33.63) 


0.12 


TTG 


6 (42.86) 


12(30.0) 


2.31 (0.62-8.63) 


0.21 



ESR1, estrogen receptor a; CDGP, constitutional delay of growth and 
puberty; OR, odds ratio; CI, confidence interval; SNP, single nucleotide 
polymorphism; rs, reference sequence. 

*ESR1 haplotypes with SNP in intron 1 arranged in the order 5' to 3'; 
rs3778609, rsl 2665044, and rs827421 . 



Table 3. Genotypicand allelic frequencies of the ESR1 in CDGP cases and controls 



SNP genotype 


Case (%) 


Control (%) 


OR (95% CI) 


/'for trend Pfor permutation 


rs3 778609; 4502C>T 








0.08 0.07 


CC 


3 (42.86) 


10(50.0) 


1.49 (0.24-9.11) 




CT 


2 (28.57) 


9 (45.0) 


2.43 (0.59-10.0) 




TT 


2 (28.57) 


1 (5.0) 


15.49 (0.8-2999) 




C allele 


8(57.14) 


29 (72.5) 


2.31 (0.62-8.63) 




T allele 


6 (42.86) 


1 1 (27.5) 


1 .88 (0.54-6.56) 




rsl 2665044; 21 187C>T 








0.22 0.18 


CC 


3 (42.86) 


9 (45.0) 


1.49 (0.24-9.11) 




CT 


2 (28.57) 


9 (45.0) 


1 .83 (0.52-6.45) 




TT 


2 (28.57) 


2(10.0) 


5.45 (0.45-66.7) 




C allele 


8(57.14) 


27 (67.5) 


1 .94 (0.53-7.08) 




T allele 


6 (42.86) 


13(32.5) 


1.67 (0.48-5.81) 




rs827421;28437A>G 








0.01 0.03 


AA 


1 (14.29) 


7 (35.0) 


4.63 (0.42-50.82) 




AG 


3 (42.86) 


1 2 (60.0) 


6.49(1.07-39.58) 




GG 


3 (42.86) 


1 (5.0) 


1915 (1.28-287.17) 




Aallele 


5(35.71) 


26 (65.0) 


4.1 (1.1-15.33) 




G allele 


9 (64.29) 


14(35.0) 


3.6(1.01-12.79) 





ESR1, estrogen receptor a; CDGP, constitutional delay of growth and puberty; SNP, single nucleotide polymorphism; OR, odds ratio; CI, 
confidence interval; rs, reference sequence. 
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a dominant effect . Genome wide association studies are 
useful method to identify the linkage of polygenic inheritance, 
and the Early Growth Genetics Consortium in Europe found 
associations of LIN28B,MAPK3, and ADCY3-POMC with 
pubertal height growth, menarche, body mass index, and early 
puberty 2 '. Also, many candidate genes were studied for a long 
time. 

Banerjee et al. 23) studied if there is any association oi leptin or 
leptin receptor polymorphism with CDGP, but no association 
was found, leptin was thought as only a prerequisite factor for 
puberty. GnRH and GnRH receptor (GnRHR) genes were also 
strong candidates because puberty begins with GnRH pulse, and 
Lin et al. 24> found a homozygous R262Q mutation in GnRHR 
in two boys with hypogonadotropic hypogonadism from one 
family. However, Sedlmeyer et al. 25) studied and concluded 
genetic variation in GnRHl and GnRHR was not likely to be a 
modulator of pubertal timing in general population. 

Pugliese-Pires et al. 26> reported the cases of CDGP with novel 
inactivation mutations in the GHSR gene, and he thought 
abnormalities in ghrelin receptor function may influence the 
phenotype of CDGP. Kauffman et al. 27) , and Seminara et al. 28) 
studied Kissl and GPR54 as candidate genes in mice. Tommiska 
et al. 29) and Ong et al. 30) also thought and studied LIN28B as 
a related gene with obesity and earlier puberty. But, Gajdos 
et al. and Vaaralahti et al. 32) 's studies told us that the gene 
variants or defects found in hypogonadotropic hypogonadism 
are not related with CDGP. Finally, it is still a conundrum to 
deciphering the genetics of pubertal onset, timing, and CDGP. 
In 2010, Borjesson et al. suggested the low concentration of 
estrogen may be involved in pubertal growth spurt and the 
interaction of high concentration of it and its receptor (ESR1) 
can be a factor of epiphyseal plate closure in long bones. 

ESR1 is in chromosome 6 and composed of 8 exons in 
human 12 '. Futhermore, there were several studies regarding the 
associations of ESR1 polymorphisms, especially in intron 1 
region due to possibly promoter effects, with height variation 
and onset of menarche 1418 '. Therefore, our hypothesis is that 
ESR1 can be a factor that controls the tempo of growth and 
puberty and a marker of CDGP. 

In this study, at the rs827421 SNP locus in intron 1 of ESR1, 
we found the frequency of G/G genotype was increased 
in CDGP patients. We thought that this means it may be a 
contributing factor of pubertal growth and a marker of CDGP. 
This is the first study about the association of ESR1 with CDGP, 
but, there is a limitation that due to the relatively small number 
of cases and a difference of numbers between the case and 
control group. Further studies are required in a large population 
cohort. 
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